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Abstract 

We propose a new method to estimate cosmological parameters using the baryon fraction of clusters of 
galaxies for a range of redshifts. The basic assumption is that the baryon fraction of clusters is constant, 
which is a reasonable assumption when it is averaged within a Mpc scale. We find that the baryon fraction 
vs. redshift diagram can estimate the cosmological parameters, since the derived value of the baryon 
fraction from observations depends on the adopted value of cosmological parameters. We also discuss some 
points in comparing theoretical calculations with observations. 
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1. Introduction 

Clusters of galaxies are the largest virialized systems 
in the universe and their masses can be estimated by 
X-ray and optical observations. In general, since there 
is no efficient way of changing the baryon fraction aver- 
aged within a Mpc scale, the baryon fraction of a cluster 
of galaxies M b /M tot (M b and M tot are baryon and to- 
tal masses of a cluster) should take the global fraction 
in the universe: flb/^o (^6 and S^o are the baryon den- 
sity parameter and the total density parameter). White 
et al. (1993) estimated fio by identifying Mb/Mtot of 
the Coma cluster with flb/flo and by using the big bang 
nucleosynthesis value for fib, and they found that low 
ft value is favored: fl = 0.16/i~ 1 / 2 /(l + 0.19/i 3 / 2 ) 
(where Ho = lOOh km s" 1 Mpc -1 is the Hubble con- 
stant). However, in this method, the derived depends 
on the Hubble constant and studies of the big bang nu- 
cleosynthesis. Furthermore, this method does not probe 
the value of the cosmological constant A . 

In this Letter, we propose a new method to estimate 
cosmological parameters (both Qq and Ao) using the 
baryon fraction of clusters for a range of redshifts. 

2. Cosmological parameter dependence on mass 
estimation 

First of all, we briefly summarize the method to es- 
timate gas, galaxy and total masses of a cluster from 
X-ray and optical observations paying attention to their 
cosmological parameter dependence. 

Firstly, we describe the gas mass in a cluster of galaxies 
using X-ray data. For simplicity, we assume that the hot 



gas in a cluster follows a spherically symmetric isother- 
mal 0— model (e.g., Sarazin 1988): 



n e (r) = n e0 1 + 



-3/3/2 



(1) 



where n e is the electron number density, and n e0 and 
r c are the central electron density and the core radius, 
respectively. Then, using X-ray data, the gas mass 
M gas (< R) within a radius R and the bolometric lumi- 
nosity Lx(< R) derived by X-ray observation are written 
as 



M gas {<R) = 
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where p gas is the gas density and other symbols have 
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their usual meanings. Combining equations (2) and (4) 
yields the gas mass: 



thus. 
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In the above equations, we implicitly assumed that the 
X-rays are emitted by the thermal bremsstrahlung only 
(see e.g., Rybicki and Lightman 1979). Note that L Xl r c 
and R are not the quantities derived by observations di- 
rectly, but that they depend on the adopted cosmological 
parameters. That is, they are estimated as 

L X (<R) = 4ir[D L (z;n ,\ ,H Q )} 2 f x (<6), (7) 
r c = c D A (z;n o ,\o,Ho), (8) 
R = 9D A (z;n ,Xo,H ), (9) 

using directly observed quantities: the total bolometric 
flux fx(< 9) within 9, the angular core radius 6 C and 
the outer angular radius 9 (see, e.g., Peebles 1993). In 
the above equations, Dl and Da are the luminosity dis- 
tance and the angular diameter distance which depend 
on the redshift z of the cluster and cosmological parame- 
ters: fioi Ao and H n . Using the relation Dl = (1 + z) 2 Da 
and paying attention to their cosmological parameter de- 
pendences, the gas mass is written as 



M gas {< 9) cx /r 5 / 2 -D/ /2 (2;Oo,A ) 



(10) 



where Da = Ii~ 1 Da- If we estimate the gas mass fixing 
the value of f^o and Aq, it is proportional to h~ 5 / 2 , which 
is a well known relation. 

Secondly, we describe the galaxy mass in a cluster. In 
general, it is estimated as the total blue (absolute) lu- 
minosity times the mass-to-luminosity ratio (e.g. White 
et al. 1993). Note that the total blue luminosity is pro- 
portional to Dl 2 and that the mass-to-luminosity ratio 
is proportional to h (since it is derived for nearby galax- 
ies, it does not depend on the cosmological parameters 
except for the Hubble constant). Thus, the galaxy mass 
Mg a i hi a cluster is written as 



M„ a i oc h x D~a 



(11) 



If we fix the value of f2o and Ao, M ga i is proportional to 

h-K 

Finally, we estimate the total mass in a cluster of galax- 
ies using X-ray data. We assume that the intracluster gas 
is in hydrostatic equilibrium. Then, the total mass within 
the radius R is 



M tot (< R) 
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3. Method to estimate cosmological parameters 

The basic assumption of our new method to estimate 
cosmological parameters is that the baryon fraction of 
clusters of galaxies Mt,/M tot is constant. In this method 
we do not need to assume that M /M to t equals the global 
baryon fraction Q /Qq as in previous works (e.g., White 
et al. 1993). Based on this assumption, we estimate the 
cosmological parameters. 

Using the results mentioned in the previous section, we 
can write the baryon fraction (Mb = M gas + M ga i) in a 
cluster as 



M b 
M tnt 



= AD A + Bh~ 3/2 D A 



~ 3/2 



(14) 



r=R 



where A and B are some cosmological parameter- 
independent values. In this Letter, we concern with only 
cosmological parameter dependence of the baryon frac- 
tion, thus precise formula of A and B is not needed be- 
low. In the above equation, the first term comes from 
the galaxies and the second from the gas. As shown in 
the above equation, the derived baryon fraction depends 
on the cosmological parameters through the angular di- 
ameter distance. 

The basic and plausible assumption adopted in this 
Letter is that the value of Mb/M tot is constant for all 
clusters. However, the derived value of M /M to t depends 
on the adopted cosmological parameters (Qq and Ao) in 
particular at high redshift, since it depends on the an- 
gular diameter distance which depends on the adopted 
cosmological parameters. That is, if we use another cos- 
mological parameters to derive Mb/M to t, we get another 
value for it. If we adopt the correct values as flo and Ao 
to estimate mass of a cluster, we obtain the true value of 
Mb/Mtot and it must be constant for all clusters at any 
redshift. On the other hand, if we adopt incorrect val- 
ues for f2o and Ao, we obtain incorrect value of Mb/M to t 
and it varies with redshift. Thus, we can determine f2o 
and Ao from Mb/M to t as follows: Estimate Mb/M to t for 
clusters with various redshifts adopting some fixed val- 
ues as fio and Ao (hereafter Mb/M to t vs. z diagram) . If 
the derived values of Mb/M to t vary with redshift, then, 
adopted values of fig and Ao are incorrect. If the derived 
values of Mb/M to t are constant, then adopted values of 
fio and Ao are the correct ones. 

It is to be noted that derived Mb/M to t depends on 
ilo, Ao and h (sec eq. (|l4|) ) . However, the Hubble constant 
dependence cannot be discriminated unless we know the 
true value of B in this method. Although uncertainty of 
the Hubble constant shifts the absolute value of derived 
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Fig. 1. The expected artificial redshift evolution of the 
derived M ga s/Mtot- The values of M ga s/Mtot are 
estimated fixing (f2o,Ao) = (0,0) when the correct 
values are (1,0) (solid), (0,0) (dotted) and (0.1,0.9) 
(dotted). Normalization of M ga s/Mtot is arbitrary. 



Mb/Mtot, it does not change their behavior of the red- 
shift evolution. Thus, we can determine flo and Ao with- 
out knowing the value of h. Furthermore, in the above 
discussion, we took some assumptions for simplicity, but 
the final result (eq.([l4|)) is correct even if we relax some 
assumptions, for example, the gas is isothermal. 

In the following, we study M gas /Mt t instead of 
Mb/Mtot, for simplicity. This is justified when the 
galaxy mass is negligible compared with the gas mass 
(Mgai/Mgas is about 20ft 3 / 2 % in the Coma cluster: 
White et al. 1993) and/or the value of M gas /M to t is con- 
stant for all clusters. In this case, we can use more sim- 
ple way to estimate Slo and Ao: We estimate Mg as /M to t 
using fixed values of Qo and Ao- We can predict the 
artificial redshift evolution of the derived M gas /M to t if 
the correct values of Slo and Ao are different from the 
adopted values. Thus, we can also determine f2o and Ao 
from derived M gas j M to t (using fixed values of f^o and Ao) 
comparing with expected artificial z evolutions. In figure 
1, we show the expected artificial redshift evolution of 
M gas /Mt t fixing (Q , X ) = (0,0) when the correct val- 
ues of (fi , A ) are (1, 0), (0, 0) and (0.1, 0.9). Of course, 
if the correct values of (f2 ,Ao) are (0,0), the derived 
values of M gas /M to t are constant. On the other hand, 
if (fi 0s A ) = (1,0)((0.1,0.9)), then derived M gas /M t ot 
increases (decreases) with redshift. If we can determine 
Mgas/Mtot within ~ ±10% for a cluster at z ~ 0.4 be- 
sides for nearby clusters, we can discriminate three mod- 
els. 



4. Discussion and Conclusions 

As mentioned above, the M gas /M to t (or Mb/Mtot) vs. 
z diagram is a potentially useful method to estimate cos- 
mological parameters. However, to estimate them by this 
method in practice, we need to check some uncertainties. 
Firstly, the value of M sas /M tot depends on the outer ra- 
dius of the analysis except for (3 — 2/3. Thus, we must 
know how large a radius is required to derive M gas /M to t- 
Secondly, we must know whether M sas /M tot (Mb/Mtot) 
is constant or not, in fact. Thirdly, we must know how 
large a scatter there is. To check these uncertainties, we 
need first study M gas /M to t (Mb/Mtot) for nearby clusters 
which are not affected by adopted cosmological parame- 
ters. Unfortunately at the present, we do not have large 
samples enough to study this, but in the near future we 
will be able to have such samples. 

At the present, we do not have large samples of 
M g as/Mtot data enough to study its cosmological param- 
eter dependences, in particular at high redshift. Thus, 
here, we use inhomogeneous samples which are collected 
from some literatures, to study the possibility of the 
method mentioned above. It is to be noted that analysis 
described below using inhomogeneous samples is not con- 
clusive, but illustrative at best. The data are summarized 
in table 1. In table 1, we assumed flo — 0, Ao = and 
h = 1. When the original data are derived using differ- 
ent cosmological parameters, we translate them using the 
above relation (eq.(|l4|)) to the values for (fio> Ao) = (0, 0) 
and h = 1. As mentioned above, this data set is inho- 
mogeneous: For example, outer radius of the analysis 
(see table 1) and method of mass estimation are differ- 
ent for each cluster and literature. In figure 2, we plot 
Mgas/Mtot as a function of redshift with expected arti- 
ficial z evolution curves. In figure 2, the values of data 
are normalized by the value of the Coma (which is the 
leftmost one with error bar). Since, here, we are inter- 
ested in only their redshift evolution and not their ab- 
solute values, the normalization is arbitrary. As is seen, 
there is a large scatter. A part of scatter comes from 
inhomogeneity of the data set. We think large homoge- 
neous data set will conclude whether this scatter is real 
or not and how large a scatter there is. If the scatter 
is real and large, we need some bayon segregation pro- 
cesses for cluster formation and evolution, for example, 
energy input by galactic winds. Since we think, in gen- 
eral, there is no efficient process of changing the baryon 
fraction averaged within a Mpc scale, we pay attention 
to data for which their outer radius is greater than IMpc 
(Oo = 0, Ao = 0, and h = 1) and they are shown as large 
crosses. The result suggests that high-Ao model seems 
to be disfavored. However, any conclusion must be left 
until large homogeneous data set is available. 

In summary, we proposed a new method to estimate 
cosmological parameters using the baryon fraction of 
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Table 1. Data for the baryon fraction. 



cluster 



redshift 



-(<#)(%) flO^Mpc) ref. 



Coma 


0.0232 


4.95 


1.5 




1 


A2199 


0.0299 


3.2 


1.2 




3 


A85 


0.0521 


6.65 


0.71 




2 


A3266 


0.0545 


5.55 


0.71 




2 


A2319 


0.0559 


4.31 


0.70 




2 


A2256 


0.0581 


5.4 


1.2 




3 


A 1795 


0.0621 


7.28 


0.71 




2 


A644 


0.0704 


3.71 


0.60 




2 


A401 


0.0748 


5.0 


1.2 




3 


A2029 


0.0768 


5.7 


1.2 




3 


A1650 


0.0840 


4.17 


0.55 




2 


A478 


0.0881 


9.05 


0.98 




2 


A2142 


0.0899 


4.99 


0.97 




2 


A3186 


0.1270 


6.54 


0.75 




2 


A1413 


0.1427 


4.07 


0.86 




2 


A545 


0.1530 


6.40 


0.91 




2 


A2009 


0.1530 


4.81 


0.65 




2 


A3888 


0.1680 


4.91 


0.56 




2 


A2218 


0.175 


6.7 


0.4 




5 


A 1689 


0.1810 


4.84 


0.74 




2 


A665 


0.1816 


7.00 


1.2 




2 


A1763 


0.1870 


7.00 


0.91 




2 


A2163 


0.2030 


5.94 


1.1 




2 


A1246 


0.216 


7.0 


1.3 




6 


C10016+16 


0.5545 


8.44 


1.7 
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The values of 


Mgas/Mto 


and R are estimated for flo 


= 0,A = 


and H = 100km s" 1 


Mpc . The last column refer to 


1. 


White et al. 


(1993). 2. 


White and Fabian 


(1995). 


3. Buote 


and Canizares (1996). L 


. Neumann and Bohringer 


(1996) 


5. 


Squires et al. 


(1996). 6. 


Kikuchi (1996). 









clusters of galaxies for a range of redshifts. We found that 
the Mb/Mtot vs. z diagram provides information about 
the values of cosmological parameters if we can deter- 
mine the baryon fraction within ~ ±10% for a cluster at 
z ~ 0.4. Furthermore, if the baryon fraction of a cluster 
is identified with the global baryon fraction as in many 
previous works, we can estimate fib as fin x (Mb/M to t) 
using fin value derived from the Mb/Mtot vs. z diagram. 
This is a new fif, estimation method without using big 
bang nucleosynthesis studies. 




. 2. Same as Fig.l. Crosses are observational data 
(see table 1) which are normalized by the value 
of the Coma (which is the leftmost one with er- 
ror bar). Large crosses indicate clusters with R > 
lh Mpc and small crosses indicate clusters with 
R < l/i _1 Mpc. We assume f2n = and An = 0. 
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